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ABSTRACT 

The coalescence of compact objects are a promising astrophysical sources of gravitational wave 
(GW) signals. The ejection of r-process material from such mergers may lead to a radioactively- 
powered electromagnetic counterpart which, if discovered, would enhance the science return of a GW 
detection. As very little is known about the optical properties of heavy r-process elements, previous 
light curve models have adopted opacities similar to those of iron group elements. Here we report 
that the presence of heavier elements, particularly the lanthanides, increase the eject a opacity by 
several orders of magnitude. We include these higher opacities in time dependent, mult i- wavelength 
radiative transport calculations to predict the broadband light curves of one-dimensional models over 
a range of parameters (ejecta masses ~ 10~^ — lO^^M© and velocities ^ 0.1 — 0.3 c). We find that 
the higher opacities lead to much longer duration light curves which can last a week or more. The 
emission is shifted toward the infrared bands due to strong optical line blanketing, and the colors at 
later times are representative of a blackbody near the recombination temperature of the lanthanides 
(T ~ 2500 K). We further consider the case in which a second mass outfiow, composed of ^^Ni, is 
ejected from a disk wind, and show that the net result is a distinctive two component spectral energy 
distribution, with a bright optical peak due to ^^Ni and an infrared peak due to r-process ejecta. We 
briefiy consider the prospects for detection and identification of these transients. 
Subject headings: radiative transfer - supernovae 



1. INTRODUCTION 

A deeper understanding of compact object binary 
mergers has the potential to shed light on several 
important and unresolved questions in astrophysics. 
Binary neutron star mergers (NSMs) may be the 
central engines of short-duration ^a mma ray bursts 
(jPaczvnskil 119861 : iNaravan et al.l Il992f ) , and heavy ele- 
ment nucleosynthesis in merger eje cta undergoing decom- 
press i on fro m nuclear densities (|Lattimer fc Schramm! 
[1973 . [T976h may contribute to the pro duction of r- 
proce ss elements in the U niverse (e.g. [Eichler et ahl 
19891: lR.osswog et al.l 119981; »Freiburghau s et al.l 119991 : 
R.osswog et al.ll200 0: Rosswog 2005; Goriel v et al1l2005D . 
With the LIGO (jAbramovici et al.. .1992. ) and VIRGO 
experiments preparing to upgrade to advanced status, 
compact object mergers provide the most likely source 
of kHz gravitational waves (GWs). The science returns 
from detected GW signals can be enhanced by the iden- 
tification of a coincident electro magnetic sign ature (e.g. 
ISchutz 1986, 2002; Stubbs 20081: iBloomlFaLl 12009). 

iMetzger fc Bergen (;2012i) consider several possible EM 
transients associated with compact object mergers - 
short-duration gamma ray bursts, orphan radio and opti- 
cal afterglows, and optical "kilonovae" light curves - and 
conclude that the last of these offers perhaps the best 
opportunity for localizing GW signals. The kilonovae 
are optical transients powered by the radioactive decay 
of material ejected in the merger. Simulations suggest 
that somewhere in the range of 10~^ — 10~^ Mq may 
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be ejected either by tidal stripping during the merger it- 
self, or later by a disk wind driven in the evolution of 
a post-merger remnant. The tidally ejected material is 
very neutron rich and should form heavy elements via 
the r-process, while the subsequent outfiows are proba- 
bly less neutron rich and nucleosynthesis may not extend 
much beyond the iron peak. In either case, the bright- 
ness of the resulting light curves will be set by the mass of 
radioisotopes ejected, while the duration depends on the 
effective diffusion time through the ejecta, and hence on 
the mass, velocities, and opacities. While the nucleosyn- 
thesis and radioactiv e heating has been modeled with a 
fair degree of d etail ( Metzger et al.l 120101 : iRoberts e t al.l 
1201 ll : iGorielv et al. 20 lU). the unknown opacities of r- 
process elements have posed a major challenge to pre- 
d icting kilonova li g ht cur ves and spectra. 

iLi fc Paczvhskil (jl998r ) were the first to study the 
r-process transients from NSM. Their simplified one- 
zone model assumed spherical symmetry and a black- 
body spectrum, and absorbed much of the physical 
complexity of the system into a set of input parame- 
ters. Assuming opacities on the order of the electron 
scattering opacity, they predicted that the light curves 
would peak in the optical/ultraviolet and rise to peak 
bolo metric luminosities of 10^^ — 10^^ ergs s~^ on a 
tim escale of a day . Mo re recent, mor e detailed, mod- 
els (jMetzger et al.l 120101 : IRoberts et all 1201 ID have used 
nuclear reaction networks to determine the radioactive 
heating rate and Monte Carlo radiation transport to cal- 
culate merger light curves, under the assumption that 
the opacities of heavy r-process elements could be ap- 
proximated by those of iron. These newer models found 
qualitatively similar results, with peak bolometric lumi- 
nosities ~ 10^^ — 10^^-^ ergs s~^ and rise times around 
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Fig. 1. — Wavelength dependent expansion opacities for ejecta 
with p = 10-13 g cm-3, T = 5000 K, and texp = 1 day. The opac- 
ity of iron is calculated using both the VALD and Autostructure 
linelists to demonstrate the reliability of the latter approach. The 
r-process opacity calculated using Autostructure data for Nd is 
in fairly good agreement with that using the VALD linelist (which 
only includes extensive line data for a few heavy elements). The 
boosted r-process opacity takes into account the diversity of species 
in an r-process mixture by assuming that all lanthanides have an 
opacity comparable to Nd. 

one day, with the colors rapidly reddening post-peak. 

In this paper, we show that using more realistic opaci- 
ties of r-process material has a dramatic effect on the pre- 
dicted kilonova light curves. VV^e use improved estimates 
of the atomic data of heavy e lements derived fro m ab ini- 
tio atomic structure models (|Kasen et al.ll2Q13l hereafter 
K13). The r-process opacities we find are orders of mag- 
nitudes higher than the those of iron group elements; as 
a consequence, we predict light curves that are longer, 
dimmer, and redder than previously thought. Rather 
than peaking sharply at t 1 day, the duration of the 
bolometric light curves can last ~ 1 week. The spectral 
energy distribution (SED) is highly suppressed in the op- 
tical, with the bulk of the energy emitted in the infrared. 
Such findings can inform observational searches for an 
EM counterpart to a GW trigger by clarifying the tran- 
sient timescales, the bands in which EM emission will be 
strongest (or have the most distinct signature), and the 
distances out to which we might expect a successful EM 
detection. 

2. OPACITY OF R-PROCESS EJECTA 

Supernova calculations suggest that for complex ions 
(e.g., the iron group) bound-bound transitions domi- 
nate other forms of opacity, such as electron scatter- 
ing, f ree-free, and photoionization (e.g. lPinto fc EastmanI 
12000^. Literally millions of lines, Doppler-broadened by 
the remnant's differential velocities, will contribute to a 
psuedo-continuum bound-bound opacity. Photons trav- 



eling through the ejecta are continually Doppler-shifted 
with respect to the comoving frame, and come into res- 
onance with multiple transitions one by one. The ve- 
locity gradient of t he remnant thus enhances the ef- 
fective line opacity (|Karp et al.l 119771 ). We account for 
this effect using the expansi on opacity formalism intro- 
duced by iKarp et al.l (ITqTTI ) and further developed by 
lEastman fc Pint J ( 19931 ), where the extinction coefficient 
is given by 

This formula represents an average over discrete wave- 
length bins, where texp is the time since mass ejection, 
Ac is the central wavelength of the bin, AAc is the bin 
width, Ti is the Sobolev optical depth of a line (|Sobolevl 
ll96Qf ). and the sum runs over all lines in the bin. The 
extinction coefficient is related to the expansion opacity 
by A^exp = <^exp/p, whcrc p is the gas density. 

To calculate Sobolev line optical depths, we assume 
that the atomic level populations are set by local thermo- 
dynamic equilibrium (LTE). This approximation should 
be reasonable in the optically thick regions of ejecta, 
where the radiation field tends towards a blackbody dis- 
tribution. In applying the a Sobolev formalism, we 
make two further assumptions: first, that the intrinsic 
(Doppler) width of lines is small compared to the ve- 
locity scale over which the ejecta properties vary, and 
second, that the intrinsic profiles of strong lines do not 
overlap with other lines. While the first condition is eas- 
ily satisfied in rapidly expanding NSM ejecta, the second 
may not be (see K13), and a non-Sobolev treatment may 
ultimately be necessary for a fully rigorous treatment of 
the radiation transport. 

The expansion opacity takes a simplified form in atmo- 
spheres where most lines are extremely optically thick 
(r ^ 1). As increases, the dependence of agxp on 
optical depth is eliminated (1 — e~^^ 1), and the ex- 
pression for expansion opacity simplifies to a sum of op- 
tically thick lines. The dependence on density and other 
determinants of optical depth are concomitantly reduced. 
Under these conditions, the number of distinct optically 
thick lines in each bin becomes the most important pre- 
dictor of ejecta opacity. An exhaustive tally of lines is 
therefore essential to accurately modeling ejecta opacity. 
Unfortunately, there is relatively little line data available 
for the heavy elements {Z > 28) expected to be synthe- 
sized in NSM ejecta. We co mpiled the hue data pro- 
vided in the VALD database (|Heiter et al. 2008), which 
includes fairly extensive data for a few heavy ions (e.g. 
Cell, Celll), but very little for most others species. 

On theoretical grounds, we expect the lanthanides 
(atomic numbers Z = 58 — 72) to contribute significantly 
to ejecta opacity, due to the complicated structure of 
their valence f-shells. This argument is illustrated with a 
simple combinatorics heuristic. The number of substates 
corresponding to a given electron orbital is roughly 

C= with ^ = 2(2/ + 1), (2) 

n\{g-n)\ 

where n is the number of valence electrons and / is the an- 
gular momentum quantum number of the valence shell. 
The number of lines should scale as , and will be much 
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greater for ions with electrons in an open / (/ = 3) shell. 
Assuming lines from the two species are equally likely 
to be optically thick, we expect atmospheres containing 
lanthanides to have a much higher opacity than atmo- 
spheres of pure iron. A similar argument could be ap- 
plied to elements of the actinide series, which may also 
be produced in the merger ejecta. However, the mass 
fractions of most actinides are predicted to be low, and 
can likely be ignored. 

Initial calculations of opacities using line data from 
VALD suggested that lanthanides have a much higher 
expansion opacity than iron (see Figure [J). However, the 
limited line data available in the VALD database makes 
it difficult to predict opacities over a range of ionization 
states. Instead, we use the theoretical lanthanide line 
data of K13 derived fr om the atomic structure modeling 
code Autostructure (BadneH""2011"). This supplies ap- 
proximate radiative data for neodymium {Z = 60) and 
a few other elements over the entire wavelength range of 
interest (- 100 A - 30, 000 A). We tested the rehability 
of the Autostructure data by comparing predicted ex- 
pansion opacities of select species, including Fe and Ce, 
to those calculated using existing line data in VALD. 
Overall, the results from both data sets are found to be 
in good agreement (Figure [1]). 

Until Autostructure models for all lanthanides are 
calculated, we use an averaging scheme to predict ejecta 
opacities based on a few representative s pecies. We re- 
cast t he ejecta composition computed by [Roberts et al] 
(|201lf ) in terms of two groups of elements: iron-like 
elements (i.e., the d-block of the periodic table) and 
neodymium- like elements (i.e., the lanthanides ). The 
line data for iron is taken from iKurucz fc Belli ([1995), 
while that of Nd is derived from Autostructure calcula- 
tions. We take the mass fraction of Nd (Fe) to be equal 
to the average mass fraction of the lanthanides (d-block 
elements). The remainder of the composition is taken to 
be calcium, which serves as a low-opacity filler with an 
appropriate ionization potential. 

Since each ion species in the ejecta has a unique set of 
strong lines, and since agxp increases with the number of 
strong lines, mixtures with a greater diversity of elements 
will have higher opacity. Since both Nd and Fe both have 
intermediate complexity for their respective blocks, their 
opacities can presumably be taken as representative of 
other elements in the same region of the period table. We 
therefore assume that each lanthanide or d-block element 
in the original compositions provides an opacity equal to 
that of Nd or Fe, respectively, and arrive at a generalized 
expression for expansion opacity 

where Z runs over the representative elements (here Fe 
and Nd), Nz is the number of elements in the block rep- 
resented by Z, and pz = Xzp^ with xz the mean mass 
fraction of the representative elements. Since lanthanide 
contributions dominate the opacity, the boosting proce- 
dure effectively increases the opacity by a factor of 14, 
the number of lanthanide species. 

Figure [T] shows the expansion opacity calcu- 
lated for typical parameters of NSM ejecta {p = 
10"^^ g cm-^, T = 5 X 10^ K, and texp = 1 day). The 



values vary with temperature and density, but in general 
our calculations predict r-process opacities many orders 
of magnitude higher than those calculated for iron group 
elements. 

3. LIGHT CURVES OF i?-PROCESS TRANSIENTS 

The surprisingly high opacity of r-process material, 
discussed in ^ has important implications for the EM 
emission from NSMs. In this section, we present radia- 
tion transport calculations using our refined opacity esti- 
mates to determine the bolometric and broadband light 
curves of r-process outfiows. Our predictions diverge 
from those of earlier studies, which assumed that the 
opacities were similar to those of iron. 

3.1. Ejecta Model 

We model the NSM ejecta as a spherically symmetric 
outfiow undergoing homologous expansion. In reality, 
the ejecta may have a hi gs:hly asymmetric, "tidal tail" 
geometry (jRosswod l2005f ). 3-D transport calculations 
suggest that this asphericity makes the emission moder- 
ately anisotropic, but doe s not qualitatively c hange the 
shape of the light curves f Rob erts et aI][20TTh . We de- 
scribe the density of the e jecta using the broken po wer 
law profile introduced by IChevalier fc Sokerl ()1989f ). in 
which density decreases as in the inner layers of the 
atmosphere and as (with n > S) in the outer layers. 
The shift from n to S occurs at the transition velocity 

vt = 7.1 X lO^C. (^5i/M)'/' cm s-\ (4) 

where £^51 is the explosion energy E/10^^ ergs, M is the 
ejecta mass Mej in units of M©, and is a numerical 
constant. For v < Vt^ the density is given by 

and an analogous expression describes the outer layers. 
The constants and (p satisfy the requirement that 
the density profile integrates to the specified mass and 
energy. We also tried using an exponentially decreas- 
ing density profile, but found that the light curves were 
mostly insensitive to the details of the density structure. 
The results presented here were generated using the bro- 
ken power law profile with (n^S) = (1, 10). 

The mass and velocities of the matter ejected dur- 
ing compact object mergers is still uncertain. Simula- 
tions of NSMs (e.g lR.osswog et al.l 119991: lOechslin et jp 
120071 : iGorielv et al.ll2011l : iHotokezaka et aIll2013D a nd of 
neutron s t ar-black hole m ergers fe.g. iJanka et al.l 119991 : 
iLed I2OOII : iRosswogI l2005f ) have found that the amount 
of material ejected is sensitive to binary type, mass 
ratio, nuclear equation of state, and the treatment of 
gravity (Newtonian v. general relativistic). However, 
taken in aggregate, simulations suggest NSMs may eject 
~ 10~^ Mq of material with velocities on the order of 
P = v/c ~ 0.1. We adopt these parameters as our fidu- 
cial model, but vary the mass and velocity scales over 
the range M = 10"^ - 1O-^M0, and /3char = 0.1 - 0.3, 
where we define the characteristic velocity by 

E = ^Mej (/3charC)' . (6) 
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TABLE 1 

Peak magnitudes and bolometric light curve properties of radioactive transients 



M ■ ( Mr-. ^ 


p 




T ^ 

-t^p,bol 


, b 

''P,bol 


iviB 




Mi 


iviH 




0.1 


r-proc 


1.2 X 10 


0.65 


-9.3 


-11.2 


-12.0 


-14.0 


10~^ 


0.2 


r-proc 


1.6 X 10 


0.75 


-9.1 


-11.3 


-12.6 


-14.1 


10~^ 


0.3 


r-proc 


3.5 X 10^° 


0.15 


-9.3 


-12.4 


-13.6 


-14.7 


10"^ 


0.1 


r-proc 


5.2 X 10^° 


4.25 


-11.2 


-12.7 


-13.3 


-15.8 


10-^ 


0.2 


r-proc 


8.5 X 10^° 


1.85 


-11.4 


-12.8 


-14.0 


-15.9 


10"^ 


0.3 


r-proc 


1.7 X 10^^ 


0.25 


-11.8 


-14.1 


-15.2 


-16.5 


10-^ 


0.1 


r-proc 


2.4 X 10^^ 


8.65 


-12.9 


-14.1 


-14.6 


-17.5 


10-^ 


0.2 


r-proc 


4.1 X 10^^ 


4.35 


-13.4 


-14.3 


-15.4 


-17.5 


10-^ 


0.3 


r-proc 


7.2 X 10^^ 


0.25 


-14.0 


-15.5 


-16.7 


-18.1 


10"^ 


0.1 




3.5 X 10^° 


0.25 


-13.0 


-12.6 


-12.6 


-11.8 


10"^ 


0.1 


56^1 


3.7 X 10^^ 


0.75 


-15.4 


-14.9 


-14.5 


-14.8 



^Peak bolometric luminosity, in ergs s ^ 

^Rise time to bolometric light curve peak, in days 

We assume the eject a to be homogenous, and composed 
of either pure r-process material or pure ^^Ni. For mod- 
els with r-process elements, all zones are as sumed to have 
the sa me radioactive decay rate given by [Roberts et alJ 
(|2QTiI) . Of the decay energy, 10% is taken to be in fis- 
sion fragments, and 90% in beta decays. Of the beta 
decay energy, 75% is assumed to be lost to neutrinos, 
with the remaining 25% split equally between leptons 
and gamma-rays. Leptons and fission fragments are as- 
sumed to be thermalized locally, while we use a radiation 
transport scheme to follow the propagation and absorp- 
tion of gamma rays. This approximate apportioning of 
the dec ay energy is based on the physical considerations 
given m iMetzger et all (j2010[ ). 

3.2. Light Curves 

We generate synthetic observables of our ejecta mod- 
els using the time-depen dent multi-wavelen gth radiation 
transport code Sedona (|Kasen et aLll2006l ). Beginning 
at an initial time of 0.1 days after mass ejection, the 
code follows the temperature and density evolution of 
the expanding ejecta, taking into account radiative and 
radioactive heating as well as cooling by expansion. The 
wavelength dependent r-process and iron opacities are 
calculated in each zo ne using the Autostructure and 
iKurucz fc Belli (|1995f ) line data, respectively. The ion- 
ization and excitation state of the gas are assumed to be 
set by LTE, and lines are taken to be completely absorb- 
ing. Sedona synthesizes the emergent spectral time se- 
ries, from which we construct bolometric and broadband 
light curves. Table [T] summarizes the EM properties of 
the models we investigated. 

In Figure O we plot the bolometric light curves of r- 
process transients for a subset of the parameter grid. For 
comparison, we also show a light curve calculated using 
iron-like opacities. Introducing more realistic r-process 
opacities clearly changes the picture significantly - the 
light curves have slower rise times and broader, dimmer 
peaks than those calculated with iron-like opacities. No- 
tably, the r-process transients do not in general have the 

1 day duration thought to be a distinguishing feature 
of these events, and instead may last a week or longer. 

That higher opacities give broader, dimmer light 
curves is expected theoretically. The duration of radioac- 
tively powered bolometric light curves is set by the effec- 
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Fig. 2. — Synthetic bolometric light curves of r-process tran- 
sients with different ejecta masses (top panel) and velocities (bot- 
tom panel). For comparison, we also show the light curve of a 
model with fiducial ejecta parameters (M = lO""^ M©, ^char = 
0.1) but calculated assuming iron-like opacities (black lines). The 
higher opacities of r-process ejecta lead to significantly broader 
light curves. The models with higher ejecta velocities correspond 
to shorter rise times and steeper declines, while those with higher 
masses have greater luminosities and longer durations. 
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Fig. 3. — Synthetic broadband light curves calculated for the 
fiducial ejecta model calculated using iron-like opacities (left) and 
r-process opacities (right). The effect of r-process opacities is to 
suppress the optical emission and shift the radiation toward redder 
bands, in particular the infrared J, H and K bands. 

tive d iffusion time through the ejected material (jArnettI 



where v is a characteristic ejecta velocity and k an ap- 
propriately wavelength- averaged opacity Because the r- 
process opacities are 10 — 100 times larger than those of 
iron, the diffusion time is significantly lengthened. The 
longer diffusion time also leads to a dimmer luminosity 
at peak, as a greater fraction of the radioactive energy 
is lost due to expansion before it can be radiated. The 
models roughly obey Arnett's law - i.e., the emergent lu- 
minosity and instantaneous r adioactive energy in put are 
approximately equal at peak (|Arnettlll979l 1 19821 ). 

The r-process light curves vary with the ejecta prop- 
erties in predictable ways. Higher mass ejections give 
a greater luminosity and longer duration, due to their 
larger radioactive mass. Ejecta models with higher ki- 
netic energies have shorter rise times, reach greater peak 
luminosities, and decline more rapidly than their lower 
energy analogues. Over the reasonable range of ejecta 
parameters considered here, the light curves exhibit sig- 
nificant diversity - the peak luminosities vary by more 
than an order of magnitude, and the durations range for 
< 1 day to as long as two weeks. 

The broadband magnitudes of the models also differ 
significantly from previous expectations. Figure [3] shows 
that, compared to a model that uses iron-like opacities, 
r-process transients output much more energy in red and 
infrared bands, with a strong suppression of the optical 
emission. We find bright, broad peaks in the J, H, and 
K bands, while the U, B, and V bands are heavily line 
blanketed and decline sharply at early times. Figure [H 
shows that the colors of the fiducial r-process model red- 
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Fig. 4. — Select colors for the fiducial model calculated using 
r-process (solid lines) and iron-like (dashed lines) opacities. The 
colors of r-process transients initially redden rapidly, then reach 
a phase of constant color characterized by a blackbody at T ?^ 
2500 K. 

den rapidly over the first day or two, and afterwards 
become remarkably constant, with the SED peaking in 
the infrared at around ~ 1 /im. Other than the very red 
color, the spectra at these phases resemble those of other 
high- velocity SNe, with a pseudo blackbody continuum 
and broad (~ 200 A) spectral features (see K13). 

The behavior of the broadband light curves can be un- 
derstood by examining the photospheric properties of r- 
process transients, since the observed SED roughly corre- 
sponds to a blackbody at the photospheric temperature 
and radius. In Figure [3 we plot the velocity and temper- 
ature evolution of the photosphere, the surface defined 

by 

T(rphot) = - / 'Ep{r)p{r)dr = 1, (8) 

'J oo 

where Hpir) is the Planck mean opacity, computed from 
our r-process line data. 

In the initial phases, the photospheric velocity and 
temperature decline steadily, refiecting the decrease in 
density and cooling of the ejecta due to expansion. At 
around ~ 2 days, however, the photospheric temperature 
stabilizes at T :^ 2500 K. Since this is close to the first 
ionization temperature of the lanthanides, this plateau 
probably refiects the sharp drop in opacity that occurs 
when these elements recombine to neutral (see K13). Re- 
combination occurs in the cooler outer layers first, and a 
sharp ionization front forms in the ejecta. Photons pass 
easily through the cooler, neutral outer layers, but are 
trapped in the ionized inner regions - the photosphere 
thus forms at the ionization front. During this phase, 
the emergent colors are roughly constant in time and re- 
semble those of a blackbody at the lanthanide recombi- 
nation temperature. Over time, the recombination front 
recedes inward, reaching the center at around 14 days. 
At this point, the ejecta is nearly entirely neutral and 
transparent. 

3.3. Uncertainties in the Opacities 
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Fig. 5. — Evolution of the photospheric temperature and velocity 
in the fiducial model, using r-process opacities. The photospheric 
temperature initially declines rapidly, but at times > 2 days set- 
tles to a fixed value near the recombination temperature of the 
lanthanides. 
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Fig. 6. — Bolometric light curves for /3char = 0-1 and a range 
of eject a masses, calculated using line data with from two differ- 
ent Autostructure models of Nd, each with a somewhat different 
energy level structure. Solid curves correspond to our preferred 
Autostructure model (opt3), while dashed curves show results 
from an alternative model (opt2). The moderate differences can 
be taken as a rough estimate of the error resulting from uncertain- 
ties in the Nd structure models. 

Though our r-process opacities represent an improve- 
ment over previously available data, the Autostructure 
models of Nd are subject to uncertainties. In particular, 
the structure models rely on an ab-initio optimization, 
such that the predicted atomic level energies and line 
wavelengths generally differ from the experimental val- 
ues. To explore the effects of these uncertainties, we cal- 
culated bolometric and broadband light curves using two 
different Autostructure models of Nd {opt2 and optS 
from K13) each with a somewhat different energy level 
structure. The optS model-which we have used in the 
calculations presented thus far-reproduces the low lying 




Days 

Fig. 7. — A comparison of broadband curves of the fiducial model, 
calculated using line data with from two different Autostructure 
models of Nd, each with a somewhat different energy level struc- 
ture. Calculations using the opt2 line data (dashed lines) predicts 
higher magnitudes for the optical B-, and R-bands than those with 
the opt3 line data (solid lines) 

energy levels of Ndll quite well, while the opt2 model 
has generally higher excitation energies which are harder 
to populate under our assumption of LTE. The result is 
fewer strong lines and a lower overall opacity for the opt2 
case. 

In Figure m we plot bolometric light curves for ejecta 
models with Pchar = 0.1 and a variety of masses. The 
light curves calculated using the opt2 line data have 
somewhat sharper, more luminous peaks and swifter de- 
chnes, consistent with the expected lower opacities. The 
differences, however, are fairly modest, and the bolo- 
metric luminosity never differs by more than a factor 
of ~ 3. The effects are more noticeable in the broadband 
light curves (Figure [7j)- In particular, the R-band light 
curves are ~ 1 — 2 magnitudes brighter for the opt2 data, 
which could have important implications for detectabil- 
ity Based on comparison to experiment, we expect the 
opts data to be more reliable, however further work re- 
fining the opacities is clearly warranted. 

In addition to errors inherent in the individual struc- 
ture models, a perhaps larger uncertainty in our r- 
process opacities arrises from the fact that we represent 
all lanthanides with the radiative data of Nd. In fact, 
the lanthanides with a nearly half open f-shell (in partic- 
ular gadolinium) will be significantly more complex, with 
perhaps ~ 10 times as many lines of Nd. We therefore 
suspect that our current calculations underestimate the 
true opacity. The light curves we present here suggest 
that increasing the opacity further will lead to an even 
greater suppression of optical emission. 

3.4. A ^^Wi-powered transient 

Given that the radioactive light curves of NSMs depend 
strongly on the composition of the ejecta, it is worth 
considering whether any elements lighter than the Ian- 
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Fig. 8. — Bolometric light curves for models that include two 
components: r-process material (from tidal tails) and ^^Ni (from a 
disk wind). We plot the light curve of the fiducial r-process ejecta 
model (Mrp = 10~^ black line) along with two models of 

pure ^^Ni with different masses (Mni = 10~^ Mq and 10~^ Mq, 
green and purple solid lines). The dashed lines give the combined 
two-component light curves. 
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Fig. 9. — A comparison of select broadband light curves for a 
pure r-process transient (solid lines) and an r-process transient 
combined with a ^^Ni-powered outfiow (dashed lines). The bluer 
SED from the ^^Ni shifts the magnitudes of the bluer bands of the 
combined SED upward relative to a pure r-process model. This 
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thanides may be produced in these events. While the 
material dynamically ejected in the merger itself (the 
tidal tails) is thought to undergo robust r-process nu- 
cleosynthesis, it is plausible that a comparable amount 
of mass may subsequently be blown off in winds from an 
accretion disk surrounding the merged remnant. Though 



the physical properties of the disk winds remain uncer- 
tain, neutrino irradiation may drive the electron frac- 
tion to Ye ^ 0.4, in which case the nucleosynthesis 
may not extend pas t Z ^ 50 (iSurman et al.l [2QQ6l. 120081 : 
iMetzger et all 120081 : iDarbha et al.l l2010f ). If Ye is very 
close to 0.5, the composition will be primarily ^^Ni. In 
this case, the EM signature of a merger may be a super- 
position of a ^^Ni- and a r-process-powered transient. 

To address this possibility, we consider a simplified sce- 
nario where 10~^ — 10~^Mq of pure ^^Ni is blown off in a 
wind immediately post-merger. Consistent with our use 
of spherical symmetry thus far, we model this wind as 
a spherical outfiow, with /3char = 0.1 and the same bro- 
ken power law density profile with (n, ^) = (1,10). We 
consider the tidal tails and disk wind to be two separate, 
non-interacting components, which is perhaps not unrea- 
sonable given that the winds are likely collimated in the 
polar regions, while the tidal tails are largely confined to 
the orbital plane. Ignoring viewing angle effects, we take 
the two component light curve to simply be the superpo- 
sition of the individual ^^Ni-powered and the r-process 
powered light curves. 

Figure [8] shows the two component fight curves, for 
two different ratios of the ^^Ni wind mass (Mni) to the 
r-process tidal tail mass (Mrp). For Mni <^ Mrp, the 
primary effect of the ^^Ni wind is to raise the early-time 
luminosity, creating a very short peak at t ~ 1 day, which 
blends into the long, fiat, r-process light curve. The cu- 
mulative light curve thus appears to have a faster rise 
time and longer plateau. If Mni ^ Mrp, the ^^Ni emis- 
sion dominates the r-process emission for the first ~ 5 
days post merger, with the two components contributing 
roughly equally thereafter. The net effect is a gradually 
declining light curve, with the long r-process plateau ob- 
scured by the ^^Ni-powered light curve. 

The addition of a ^^Ni component also affects the 
SED of the transient, as shown in Figure [9] for the case 
Mni = Mrp = 1O~^M0. Given the much lower iron group 
opacities, the SED of the ^^Ni ejecta is much bluer than 
that of the r-process ejecta. The emission in the op- 
tical bands (U,B,V,R) is relatively bright and set by 
^^Ni mass, while the r-process material establishes the 
behavior in the infrared bands. Such an unusual SED 
may serve as an EM fingerprint that could improve the 
prospects for positively identifying a NSM. In particular, 
as shown in Figure [lOl the spectrum of a two compo- 
nent outfiow is, to first approximation, the superposition 
of two blackbodies - a sharply peaked bluer blackbody, 
corresponding to the ^^Ni ejecta, and a lower, redder one, 
corresponding to the r-process material. 

The aggregate light curve model we present here 
glosses over some of the more complex physical processes. 
Our model assumes spatially distinct regions of pure ^^Ni 
and pure r-process material. In reality, the nucleosyn- 
thetic yields are highly sensitive to the conditions in 
the wind, and it is possible that disk outfiows contain 
some elements heavier than ^^Ni. Contamination of the 
outfiows with even a small mass fraction of lanthanides 
(~ 10~^) can significantly increase the opacities and the 
optical line blanketing. Even if our simplified compo- 
sitions turn out to be reasonable, our model does not 
account for the geometry of the ejecta and any possible 
mixing of the wind and tidal tail components. Given 
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Fig. 10. — A combined ^^Ni and r-process spectrum at t = 7 
days, taking Mni = Mrp = 10~^Mq. The peak at blue wavelengths 
is due to the ^^Ni while the r-process material supplies the red 
and infrared emission. The best fit blackbody curves to the indi- 
vidual spectra are overplotted in dashed black lines (Tni — 5700 
K, Trp ~ 2400 K). The combined spectrum roughly resembles the 
superposition of two blackbodies at different temperatures. 

the presumably high level of asymmetry, the net (tails 
+ wind) EM output may depend heavily on orientation, 
making our simple superposition procedure valid only 
along certain lines of sight. 

4. CONCLUSION 

We have shown that the radioactive powered light 
curves associated with NSMs are greatly modified when 
more realistic values for the opacities of r-process mate- 
rial are taken into account. The r-process opacities are 
much higher than those of iron, due to both the com- 
plexity of heavy elements (in particular the lanthanides) 
and the diversity of atomic species present. Refining our 
understanding of the atomic structure of these elements 
is an important step toward a more rigorous model of 
transients from merging compact objects. 

In accordance with theoretical expectations, the ex- 
tremely high r-process opacities result in bolometric light 
curves that are broader and dimmer than those calcu- 
lated assuming iron-like opacities. Our calculations in- 
dicate that the light curves are likely to last at least a 
few days, and may endure as long as a week or two in 
certain cases. The broadband magnitudes are also signif- 
icantly impacted - we find heavy line blanketing in the 
optical and UV bands, with most of the radiation emit- 
ted in the near infrared. The colors at later times are 
fairly constant, and regulated to be similar to a black- 
body at T ^ 2500 K, the recombination temperature of 
the lanthanides. 

These findings have important, if mixed, consequences 
for the detectability of EM counterparts to NSMs. On 
the one hand, we predict dimmer bolometric luminosities 
and SEDs largely shifted into the infrared, both of which 
pose serious challenges to observational surveys at opti- 
cal wavelengths. On the other hand, the light curves are 



of longer duration, and may not require quite as a high 
cadence of observations. Perhaps more importantly, the 
uniquely high opacity of r-process ejecta provides signa- 
tures that may allow us to distinguish NSMs from other 
sorts of dim transients. In particular, the SED of r- 
process ejecta peaks in the infrared, with a color tem- 
perature set by lanthanide recombination. If the merger 
ejects two separate mass components - r-process tidal 
tails and a ^^Ni wind - the dual spectrum may be quite 
distinctive, with discernible infrared and optical compo- 
nents. 

The SEDs we predict can be used to roughly esti- 
mate the detectability, given the varying depths and 
wavelength coverage of different observing facilities (e.g., 
Nissanke e t all 120121 ). For example. Pan- ST ARRS (see 
http : / / pan-starrs.ifa.hawaii.edu) and PTE (|Law et al.l 
?Q09l ) achieve an R-band depth of Mji ~ 21 mag- 
nitudes, while LSST reaches a d epth of Mr ~ 24 
(|LSST Science Collaborations! I2009D . We find that an 
r-process transient with fiducial model parameters will 
peak at Mn = — 13, which under ideal observing condi- 
tions, would be observable to Pan-STARRS or PTE out 
to a distance of ~ 60 Mpc. This is an interesting, but 
rather small fraction of the volume probed by advanced 
LIGO/VIRGO. The case with LSST is more promising, 
with sensitivity in the R-band out to ~ 250 Mpc. Dis- 
covery of r-process ejecta in the U or B bands with any 
facility would appear to be quite difficult, given the heavy 
line blanketing at these wavelengths. 

Given that our models predict that most of the emis- 
sion is at longer wavelengths, improving detection ca- 
pabilities in the near infrared may greatly aid in future 
searches for EM counterparts. Ground based facilities 
with sensitivity in the I or Y bands (0.8 — 1.1/im) may 
benefit from these capabilities, as the r-process tran- 
sients are generally ~ 1 magnitude brighter in these 
bands than in R-band. Th e construction of space based 
facilities such as WFIR ST (|Green et al.ll2012[ ) and Euclid 
(|Amendola et al.|[2012f ) would be of particular interest. 
WFIRST is proposed to have an H-band depth of ~ 25 
magnitudes, with Euclid achieving a similar sensitivity. 
As our fiducial model is much brighter in the infrared 
{Mh — -15) than in the optical bands, such facilities 
could potentially make a detection out to a distance of 
^ 1000 Mpc, encompassing the entire LIGO/VIRGO vol- 
ume. 

Discovering the EM counterparts to NSMs would be 
made significantly easier if, in addition to r-process ele- 
ments, these events also separately eject some significant 
amount of ^^Ni or lower mass {Z < 58) radioactive iso- 
topes. Our models predict that such "lanthanide-free" 
light curves are reasonably bright in the optical bands 
{Mb ^ Mr ^ — 15) and would be within range for many 
upcoming optical transient surveys. It is plausible that 
winds from a post-merger accretion disk may produce 
such lighter element outflows, although more detailed 
simulations are needed to constrain the mass and com- 
position of the material ejected. Clearly any detection of 
a short-lived optical transient should, if possible, be im- 
mediately followed up at infrared wavelengths to look for 
a coincident r-process transient from the tidal tails. Dis- 
covery of such a two component light curve and spectrum 
would be a very strong signature of a NSM. It would also 
provide insight into the merger and post-merger physics 
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by separately constraining the mass ejected by different 
mechanisms. 

The work we have presented here is an important 
step towards improving our predictions of the radioac- 
tive transients from NSMs. However, much remains to 
be done. The opacities we have used, while more real- 
istic than previous estimates, are still subject to impor- 
tant uncertainties. In particular, we need comprehensive 
structure calculations to derive radiative data for all lan- 
thanides. In addition, more detailed simulations of the 
dynamics and nucleosynthesis of the mass ejection are 
needed to better predict the mass, composition, and ge- 
ometry of the ejecta. Of special interest are the prop- 
erties of the disk wind and any mixing of these outflows 
with the tidal tail material. Finally, 3-dimensional radia- 
tive transfer calculations will be needed to predict the 
light curves of multi-component mass ejections and to 
determine their dependence on viewing angle. Such the- 



oretical work should improve our understanding of the 
EM counterparts to gravitational wave sources, and the 
heavy elements they produce. 
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